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• Abstract: We study light-cone distribution amplitudes of heavy-light systems, such as a i3-meson. 
. By an explicit computation, we determine how two-parton distribution amplitudes mix with three- 
', parton ones at one loop: (p+ is shown to mix only into itself, whereas (/>_ mixes with the difference of 
I three-parton distribution amplitudes — ^v- We determine the corresponding anomalous dimen- 
sion and we check the gauge independence of our result by considering a general covariant gauge. 
Finally, we comment on some implications of our result for phenomenological models of these distri- 
5_i ■ bution amplitudes. 
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1. Introduction 

B physics provides outstanding opportunities to test the Standard Model. At the end of the era of 
i?-factories, a remarkable sample of processes has been measured accurately. On the theoretical side, 
the dynamics of heavy-light systems is naturally dominated by two scales : the heavy-quark mass mi, 
corresponding to hard contributions and a hadronic scale A related to soft, hadronic, contributions. It 
is therefore natural to deal with i?-mesons in the framework of the heavy-quark expansion, expanding 
observables in powers of l/rub, and to attempt to factorise hard and soft dynamics. 

This has been progressively achieved within the frameworks of QCD factorisation [1-3] and Soft- 
Collinear Effective Theory [4-8]. The hadronic input collecting soft physics are not only the form 
factors, but also the light-cone distribution amplitudes, defined generally as the matrix element of a 
non-local operator along a light-like direction. This quantity provides the amplitude of probability of 
finding partons inside a hadron with given fractions of the hadron momentum. For light hadrons, one 
can define a twist, corresponding to the power of the hard scale to which the distribution amplitude 
will contribute to a given process, which is related to the number of partons and the Lorentz structure 
of interest [9]. For instance, the leading-twist distribution amplitude of the pion can be defined as: 

{7T{pMz2)[z2,z,]^^j,q{z,)\0) = -iUp, /'dxe*(-P-2+^P-i),/.^(x) (1.1) 
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where {z2 — zi)'^ = 0, x 
the expression, reads: 



X, and the path-ordered exponential, ensuring the gauge invariance of 



[z2,zi] = Pexp 



i9s 



(1.2) 



This distribution ampUtude corresponds to picking up a pair of valence quark and anti-quark in the 
pion, carrying a fraction x and x of the pion momentum respectively. 

In the early days of QCD and the parton model, distribution amplitudes of light mesons were 
identified as key tools to analyse exclusive processes at high energy through factorisation, and they 
have been analysed extensively [10-12]. Phenomenologically, sum rules and lattice have been used to 
determine the values of their lowest moments [13-16]. Mathematically, the structure of the distribu- 
tion amplitudes was described by exploiting properties of the conformal group, providing a geomet- 
rical definition of twist compatible with the phenomenological one in the case of light mesons [9] . 

The case of heavy-light mesons has been discussed only recently [17, 18], mainly due to their 
importance in relation with i?-physics, as shown in the case of non-leptonic decays [1-3], scmileptonic 
decays [19,20], radiative decays [21-25]. It is possible to define two two-parton distribution amplitudes 
(f)+ and (j)- from the most simple non-local matrix element: 



{0\qfsiz)[z,0]{hv)amB{p)) = - 



(1.3) 



with the usual definitions of Heavy Quark Effective Theory (HQET) for the velocity v = p/Ms, the 
heavy-quark projection hy and the decay constant /b- The Fourier transforms of the distribution 
amplitudes, depending on t = v ■ z read: 



(1.4) 



It turns out that only one distribution amplitude, ^+ , enters most of the computations considered in 
the framework of factorisation (non-leptonic decays, B — > Vj, B jiu) at leading order in I/uib- 
In cases where at least one of the outgoing particles has not a light-like momentum, factorisation 
may still hold, but the formula involves the two distribution amplitudes of the B meson. This is in 
particular the case for B Vj*, which has an important potential to test the Standard Model [26,27]. 

Another interesting place where the distribution amplitudes naturally occur is light-cone sum 
rules. These sum rules allow in particular for a determination of form factors of phenomenological 
relevance such as B ^ tt, B ^ p. The purpose of sum rules is to reexpress those form factors 
in terms of an integral of the distribution amplitude of one of the external mesons with a kernel 
resulting from the expansion of a properly chosen correlator along the light-cone. Recently, sum 
rules with interesting properties have been proposed to relate these form factors with B meson 
distribution amplitudes. In this context, plays a dominant role and its modeling can improve the 
determination of the form factors [28-31]. 

One can also define distribution amplitudes beyond the two-parton level. Interestingly, one 
can exploit quark equations of motion in order to relate two-parton and three-parton distribution 
amplitudes, which are defined in the case of heavy-light mesons through: 



{0\qf3{z)[z,uz]gG^,{uz)z''[uz,0]{hM^)\B{p)) 
fB{li)M 



(1.5) 



1 + i>) {v^i^ - t-i^) {^A{t, u) - §y (i, «)) - ia^^z^'^vit, u) 



-z^XA{t,u) + ^YA{t,u) 
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A good model for the distribution amplitudes of (f)^ and must embed as many theoretical 
constraints as possible. Not much can be said for sure, apart from sum rule estimates [32,33], the 
expected behaviour at the origin ~ lu and (p-ico) ~ 1 for a; 0), and the properties of these 

quantities under renormalisation, which can be derived from perturbation theory. One must notice 
that differences between light-meson and heavy-light meson in terms of the renormalisation proper- 
ties [33] exist: the limit in one case is the chiral limit, which modifies only long-distance properties, 
whereas the heavy-quark limit affects short-distance features of the theory (the UV structure of 
HQET is qualitatively different from QCD). This explains the non-commutation of the heavy-quark 
limit with the light-cone limit (contrary to the chiral limit). In the case of heavy-quark distribution 
amplitude, there is no equivalent of Gegenbauer moments of light-meson distribution amplitudes, 
which mix only into themselves under renormalisation. 

The RGE behaviour of ^+ and ^_ has been investigated in refs. [34] and [35]. It was shown in 
particular that RGE generates a radiative tail, leading to a divergence of positive moments of these 
quantities. In particular, there is no absolute normalisation of (/)± to 1 from its zcroth moment, so 
that a naive partonic interpretation like in the pion case is not possible after renormalisation. To 
obtain the first inverse moments of (f).^ and (p-, which are relevant phenomenologically, one needs 
the knowledge of the whole distribution amplitude. But the previous behaviours and models were 
derived in the two-parton approximation, even though the equation of motions indicate the potential 
mixing with three-parton distribution amplitudes. 

The goal of this paper is to understand the renormalisation of these objects at the first nontrivial 
order of the strong coupling constant, including the contribution from three-parton distribution 
amplitudes. This task requires us to consider a three-parton external state (a quark, an antiquark 
and a gluon). In Sec. 2, we recall the behaviour of and ^_ under renormalisation as determined 
from a two-parton external state. In Sec. 3, we give the one-loop diagrams and their ultraviolet 
divergences in the case of a three-parton external state. In Sec. 4, we determine the mixing of ^+ 
and (f)- with three-parton distribution amplitudes at one loop, and we show that there is no mixing in 
the case of whereas 0_ mixes with the difference — ^v- In Sec. 4, we extend our calculation, 
performed in the Feynman gauge, to a general covariant gauge, which provides a check of the gauge 
independence of our results (more detailed results are given in App. A). In Sec. 5, we make a few 
comments before concluding. 



2. Two-parton S-meson distribution amplitudes 

We set up our framework by introducing notation and definitions, and we recall the results obtained 
on the renormalisation of the 5-meson distribution amplitudes. We define light-cone directions by 
two vectors 

n^ = n?. = n+-n_ = 2 v = {n+ + n^)/2 (2.1) 
so that an arbitrary vector can be projected as 

= (n+ • <Z)^ + (n_ ■q)^+q±,= + + (2-2) 

The computation of the renormalisation properties of the distribution amplitudes requires us to 



consider matrix elements of the relevant operators 

Of(c^) = ^/ dte}'^\{)\q{z)[zMi+'^KmH) (2.3) 

O'liu:) = ^ j dte''^'m{z)[z,Q]i_TKmH) (2.4) 

O^{io,i) = J dte'^' j due'^^\mz)[z.uz]gsG^,{uz)z^[uz,Q\TKmH) (2.5) 

with z parallel to n+, i.e. z^ = tn^^^, t = v ■ z = z_/2 and the path-ordered exponential in the n+ 
direction: 



[z,0] = Pexp 



igs / dy^Af'iy) 







(2.6) 



pi pi pa 

= l + igs I da z^A^{az) - gl l da l d^z^ z^ A^'iaz) A" {(3z) + . . . (2.7) 
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We define the different distribution amplitudes in momentum space through their Fourier transforms: 

<P±iuj) = dte"-'4>±it) F{u;,0 = j dt j dute'^'-+<^'F{t,u) (2.8) 

where F = ^y, ^A, Xa, Ya- 

The renormalisation group equation of and (p- can be determined by computing the mixing 
terms Z defined in the following expression (since we work in the chiral limit, we omit mixing terms 
between cp^ and (p- which were shown in ref. [35] to be proportional to the mass of the light-quark.) 

Of''^'"(a;,^) = duj'Z^\uj,u';fi)0^''"''%u') + J dJ d^' Z^l{oj,J ■ fi)0^'''^'\J ,£,') + .. . (2.9) 

where the ellipsis involves matrix elements of operators corresponding to a higher number of partons 
(a similar equation for operators related to a higher number of partons instead of 0± could be 
written, with operators corresponding to an arbitrary number of partons on the right-hand side). 
The behaviour under renormalisation being a short-distance property, any choice of H is allowed in 
principle. For instance, Z± were computed in refs. [34, 35] using a two-parton external state with 
on-shell quarks H = h{p)q{k), for which one obtains at leading order: 

(cj) = 6{uj - k+) vijTu 0^{uj, to') = (2.10) 

Computing NLO terms provides Z± in eq. (^]^). However convenient, this choice of external state 
prevents us from determining Z±^3 describing the mixing between 0± and O3, since vanishes 
then. 

One can still write down the RGE restricted to (j)± as 

(t>±{u),^l) = - du'j±{LJ,Lj',fl)(j)±{uj',fi) + ... (2.11) 



dlog n 







/dZ^^ {uj Cj fi^ 
doJ — ^(ii^g^' — Z±{oj,uj',fj,) - 'yFias)5{uj - lo') (2.12) 

where 'jp is due to the normalisation of the distribution amplitudes which involves the (renormalisation- 
scale dependent) HQET decay constant. 




HO LO HL 

Figure 1: Diagrams evaluated to determine the renormalisation constant Z± using a two-parton external 
state (the two letters describe the ends of the gluon line: Heavy quark, Light quark or Operator, i.e. the 
path-ordered exponential) . 



^±l/fO = ^x2x-y^ 17Va5' [<5(^-^'-/+)- 5(^-^0] (2.13) 
^±Ilo = ^x(-2)x- f ^^+^[5(a;'-a;-/+)-5(u;'-u;)] (2.14) 



Working in dimensional regularisation with d = 4 — 2e dimensions, in the Feynman gauge and 
in the MS scheme, one obtains from the three diagrams shown in fig. |l] [34,35]: 

X 2 X - / -r- { -. 

An e Jq 1+ \fi- 

1 r'd/,/_ 

4vr ^ ' eJo I. 
Z-\hl = ^><2xi r '^5{J-u:-W) ^+Ihl = (2.15) 



Only the diagram where the gluon line connects the two external legs differs for and (/>_. In 
addition, one has contributions from the wave-function renormalisation of the external legs: 

Z±Lfr-H = Z]!''='^^U{^-J) (2.16) 
^^Ur-. = = ^ X (-1) X \5[. - J) (2.17) 

so that Z± is 5{uj — u>') corrected at one loop by the sum of eqs. (|2.13| )- (|2.17| ). 
This leads to the following anomalous dimensions: 

7i" =75."-r™,^<^ (2.19) 

LO 

with 

r^ilp = 4C^ 7^'^ = -2Cf 7j,^) = -3Ci. (2.20) 

Quantities with a superscript (1) must be multiplied by j^. The anomalous dimension is of the 
Sudakov type, which is related to the fact that the operators of interest can be seen as containing 
two Wilson lines, one from the heavy quark along the v direction (representing the interaction of 
soft gluons with h^) from — oo to 0, linked with another one along the n+ direction from to z. The 
presence of a cusp between the two Wilson lines is responsible for the appearance of a Sudakov-like 
behaviour of the anomalous dimension [34,36,37]. 



3. One-loop computation for three-particle external state 



We perform the same computation, taking as an external state H = h{p)g{e,q)q{k) containing three 
partons. We compute the diagrams at one loop, using dimensional regularisation, and we pay a 




A : -Qs— [5 (uj-k+- q+) - 6 (uj - k+)] vd_rT''u 




B : -gs^S (lo - k+) vdJ^T^'u C : Qs 
V ■ q 



{k + q) 



;6 {uj-k+- q+) v^{^ + 4)4_rT''u 



Figure 2: The three leading-order contributions to the matrix element of 0± with a three-parton external 
state. 




As^ : igsiq+e^ - q^e+) vTT''u6{u; - k+)6{C - q+) 



Figure 3: Leading-order contribution to the matrix element of O^fj, with a three-parton external state. 



special attention to separating e-poles related to UV divergences and IR divergences carefully. We 
then identify the e-poles corresponding to UV divergences as part of the renormalisation function in: 

Of'*«-(^) = I d^'Z±(u;,a;';/i)Of''^^"(a;',/x) + J dio'dC Z±^3{uj,uj' , f^)0^'''"''{u^' ,e , f^) (3.1) 

Subtracting the contribution from Z± (S) will yield the mixing term Z±^^ between two and three- 
particle distribution amplitudes. 

At leading order, we obtain three different contributions, denoted A, B, C, for the matrix element 
of 0±, shown in fig. 0. The expression for C can be simplified using the following relations 



€ ■ q = "^u = u = 



(3.2) 



For O3, we have a leading-order expression indicated in fig. ^. 

One can easily spell out the diagrams for 0± at one loop by taking the above three diagrams 
ABC, and adding a new gluon line in all possible ways (our naming scheme reflects this idea). 
Diagram A yields the diagrams in fig. ^, whereas the other diagrams coming from B and C are show 
in fig. |5[ Let us notice the presence of redundant diagrams, namely {B2A) = {A12), (C24) = (A23), 
{B25) = (BU), (C25) = (C12). 

3.1 Common contributions 

Since we are interested in the renormalisation properties of the distribution amplitudes in the MS 
scheme, we quote here only the poles in s defined as d = A—2e corresponding to ultraviolet divergences 



All 



A12 



A13 



A14 




A22 



A23 



A24 



A33 




A34 



A44 



Figure 4: One-loop diagrams obtained from diagram A. 



(we discuss our integration procedure on one example explicitly in App. For the integrals going up 
to infinity, we keep the expression of the kernels before picking up the pole in e, since the integration 
may give rise to double poles in the expression of 7±,3, related to Sudakov logarithms. We do not 
give the expressions corresponding to the wave- function renormalisation of the external legs, i.e., 
(11), (22), (33). 

The following diagrams yield contributions of the same form for both distribution amplitudes : 



(^12) = -^Ca 9se+ v^^m - f 
(A14) = --l(g^e+) v4^m - 



1 

dl+ — — [5 -k+ + 1+) - 5{uj- k+)] 

oo / ;2 \ 

dL 



a. 

2Cf 
1+Q+ 

— 7 {^-k+-q+- /+) -6{uj- k+)] --^ [5 {uj-k+- 1+) -5{uj- k+)] 



{[5 {uj - k^ - q+ - /+) - 6 {uj - k+ - q^)] + [6 {uj - k+) - 6 {uj - k^ - /+)] 



(A24) = -^Ca^^ v4^m - f 



47r 



dL 



1+ 



21+ + q+ 



-^[6{uj-k++ /+) -6{uj- k+)] + — [6 {uj-k+- q+) - 6 {lo - k+)] 
^/lQ/l^ \ -J 1 f^+dl+k+-l+ 



2Cf 
Q+ 



{[6 {uj - k+ - q+) - 5 {uj - k+ - q+ + /+)] + [6 {uj - k+ + /+) - 6 {uj - k+) 



-Ca 



1 



[5 {u- k+ — q+) - S{io -k++ 1+)] 



1+ 



[6 (w — A;+ — g+ + Z+) — 5 {lo — k+ + 



(AU) = 
For diagrams of B-type, we obtain: 



(B12) 


= 


(BIS) 


= 


(BU) 


= 


(S15) 


as 1 


(534) 




(535) 




(B44) 


= 


(545) 


= B^Z±\jjo 


(555) 


= BxZh 



V ■ e 
( 

V ■ q 



For diagrams of C-type, we obtain: 

as 



47r^"2 (k + g)2 ""^^^ + ^)^±rT«n -6 (u - k^ - g+) 



(Cf - CaII) ,,^ + m^^T'^u U iLO-k+- q+) 



(C34) = C0 Z±|^o 
(C35) = Z±|^^ 
(C44) = 
(C55) = C xZq 



3.2 



The remaining diagrams yield different contributions for 0+ and <^_. For ^+ we have the following 
contributions: 



(A13+) = 

igse+) 



(A23+) 

(523+) 
(C13+) 

(C14+) 
(C45+) 



-Ca- 



1 



1 



dL 



{k+ -q+- 2/+; 



JO 



ik+ + q+)il+ + q+) 



^ U-q+ Q+ Jo 

{S{ui — k+ - q+) — S{iu — k+ + /+)} 








n 1 
— {Ca- 2CF)gs-v4+rT^ue+ 
47r £ 



Jo K Jk+ 



^=++9+ dl+ k++q+- 



q+ 



k+ + q+ 
c®z+\^o 



{5{uj - fc+ - q+) - 5{uj -k+-q+ + /+)} 



3.3 



For (^_, the remaining diagrams yield the following contributions: 

(T 1 

(^13_) = ^{Ca- 2Cf) gse+ vi^_TT''u - 



X / dL 



J 

Jo 



[S {uj -k+ — q+ + 1+) - S{uj -k+ + /+)] 



X i e+ v7^_rT''u 



r+ di+ 

Jo k+ 



dl+ 
Q+ 



U-q+-2k+}^l + 
l+ + k+- 2k 



k+ - 1+ 

k++q+ 

l+ + q+ 



k+ + . 



1 



dl+ 1+ + q+ 
k+ k+ + q+ 
^+ dl. 



+ 







dl+ {l++q+){k+ + 2q+y 



dl. 



J —i 



q+ik+ + q+) 



ik+ + q+){q+ + l+) 



J-q+ Q+ 

[6 {uj -k+ - q+) - (5 (w - /£+ + /+)] 



(S23_) = 



(C13_ 



{CU. 



(C45_ 



X 

as 



dL 



1+ 



q+{q+ + k^ 



-I 



^ {CA-2CF)gse+vi_TT''u 
An £ 



— 7 ; — r 

K{<l+ + k+) 

1 







dh 



1 



^{CA-2CF)gs- 
47r e 



dL 



k+{q+ + 



^ (a; — /c+ + Z-i 



d{uj-q+- 1+) 







k. 



kj- - L 



1 



Jo 



k+ + q+ 

^+ dl+ k+ - 1+ 
k+ k+ + q+ 



g+ 9+ V9+ k+ + q+J_ 



^ L^O ''^+ 5+ 



{k+-l+){k+ + q+) 



[S {uj-k+- q+) -5{u-q+- /+)] 



/ 

^0 



L 



1+ {k++q+Y 



[5 (uj-k+-q+ + 1+) -5{u-k+- 



(C45_) contains an additional contribution compared to the two-parton case, because the de- 
termination of Z-\j^Q relied on the fact that the light-quark coming out of the vertex was on shell, 
which is not the case for (C45) in the three-parton case. 



4. One-loop mixing of 0± with 3-parton distribution amplitudes 



4.1 0+ 

Having calculated the divergent part of all possible diagrams the renormalisation matrix can be 
determined in a similar manner to ref. [17]. We write: 

Z±(a;, J- f,) = 5{u, - J) + ^zl^ /x) ^±,3(0;, J, i';ii) = i^^^' ^i''. t^). (4.1) 

with being proportional to C^?. One can schematically write for the matrix element of the bare 
operator up to one loop: 

{Q\0±{oj)\Hf'''^ = z]l'^Zy'^zl''^Zg[A + B + Cf"'^ (4.2) 
+ [BU + B35 + Bib + C34 + C35 + C45] + [512 + 515 + C12 + C15 + 555 + C55] 
+ [y412 + ^13 + ^14 + A23 + A2A + AM + 523 + C13 + C14] 

= [A + B + C]^^"(^) (4.3) 



where the renormalisation constants Z/^, Zg, Z3 and Zg come from the heavy-quark, light-quark and 
gluon external legs and the coupling constant respectively in the leading order contribution. Since 
the matrix element of the renormalized operator 0±{uj; /i) must stay finite for e — > and since we 
know z"^ {uj,uj'; fi), we can determine z^^i^i'^'^''! 1^) from the poles of the diagrams listed in (^.21). 

In the case of 5 and C, the diagrams (534), (535), (545) and (C34), (C35), (C45) together 
with the fermion wave function renormalisation Z^^"^ and Z^"^ , given in eqs. ( p.l6| )-( p.l7| ), add up 
as indicated above to 5 (g) and C ® respectively. The combination of the renormalisation 
constant for the coupling constant and the gluon field tensor is: 

= 1 - ^ (4.4) 

so that its contribution multiplied by 5 and C cancels the C^-part of (512) + (515) and (C12) + 
(C15), whereas the CiT'-part of the same diagrams is cancelled by (555) and (C55), as expected from 
general arguments on the renormalisation of the quark-gluon vertex. 

The remaining diagrams in eq. (4.2) must be added, and one has to subtract A® z^^ to extract 
the three-parton contribution, which amounts to subtracting: 



(4.5) 



[5 {uj — — — /+) — 5 {to — — q^) — 6 {to — k^ — /+) + 5 {lo — k^)] 



A® Z±\^o = 2^^r 0^ X vi^TT^'u 

dl+ 1+ + k+ 



(4.6) 



1+ -Q+q+ + 



[5 (to — k^ — Z+) — 6 {uj — k^ — g+)] 







d/+ /+ + k^ 
h -Q+Q+ + 
dl+ /+ + fc-i 
U k^ 



+ / , " " ; [5iu;-k+-l+)-6{u;-k+ 
r>0 _ ^ 

[5 {uj — k^ — /_|_) — 5 {uj — k^) 



\HL 



2^-^- X {-gs)— v4_m 



9+ 



dU 



1 



+ 







dl 



k+ k+{k++q+) 



q+ + k+ 

6 {uj — — 1+) 



6 {uj — k+ — /+) 



(4.7) 



(4.8) 



Let us focus on 0+ in the remaining part of this section. For the part proportional to Ci?, the 
Sudakov-hke contribution of (^414) matches that of Z^q and more generally, an explicit computation 
shows that the diagrams in the fourth bracket of eq. ([4.2| ) add up exactly to the contribution from 
z^^ (which is proportional to Cp). 

For the Ca part, the contribution from (^414) may seem surprising at first glance, since it seems 
to involve another Sudakov-like integral, with a double pole in 1/e generated by the integration of 
up to infinity. But let us split the first integral in the Ca term of (^14) in two intervals, from 
to q+ and from g+ to 00, perform a change of variable ^+ ^ /+ + (/+ and add the second integral, we 
obtain: 









dl+ 


1 










1 /■'?+ 




-Jo 




dl+ 









2e 



l+ + q+ 



2e 



1 



l+ + q+ 



[6 {uj-k+-q+- /+) -6iuj- k+)] (4.9) 



2e 



[6 (u - k+ - 1+) - 6 {u - k+)] 



One can see that the first integral yields no pole in 1/e, and the second integral provides a simple 
pole. Summing up all the contributions proportional to Ca (including that of the renormalisation 
constants from the gluon field tensor and the strong coupling constant), one observes that they cancel 
exactly. 

In summary, the determination of the renormalisation properties of 0+ at one loop with a three- 
parton external state yields a Cp term equal to the self-mixing obtained from the consideration of 
two-parton external state, and no Ca term. This shows that 4>^ mixes only with itself, and not with 
three-parton distribution amplitudes, up to one loop. 

The fact that (p^ occurs in most of the factorisation analyses for B-meson decays suggests that 
it holds a special status with respect to other S-meson distribution amplitudes, which is somehow 
confirmed by our finding of an absence of mixing. For light mesons, conformal symmetry would 
naturally explain the absence of mixing between distribution with different parton numbers and thus 
different twists [9]. In the heavy-quark case, conformal symmetry cannot be invoked anymore [33], 
but our result may be the hint of another symmetry singling out 0+ with respect to other distribution 
amplitudes and explaining that 7+^3 = 0. 



4.2 (/)_ 

We consider now (j)^. One can follow the same argument as before, with a similar pattern of cancel- 
lation for the diagrams yielding the same results for and (/>_. In particular, one recovers easily the 
contribution proportional to Z_, i.e. the contribution from self-mixing derived from the two-particle 
case. 

But the diagrams do not cancel completely, and there remains a genuine three-particle contribu- 



tion: 



(0|O_(a;)|iJ)^"'-^ = [^ + 5 + Cr"(Ai) + ^ J (L;'z^2\^,^'; f^)[A + B + C]icj') 



1 n 1 

X |(Ca-2Cf) 



qI Jk. 



1 „ / 1 1 \ 1 

^ [h'Kj^ WT^f Jo 

1 r 1 rH+g+ I rq+ "i ^ 



X {5{uj -k+-q+ + 1+) - d{uj -k+- q+)} , 
trate tl 



from which 7pi^^ can be extracted. If we separate the Dirac structure 



(4.10) 
(4.11) 



we obtain 
.(1) 



% 

'2e 



CO 



- Ca 

We defined the +-distribution as: 



(Ca - 2Cf) ,2 / , tf 

OJ' + — u 

e{uj' + -to) 1 



/2 



{to' + eo' 

{Q{lo - u}') - Q{lo -J - CO) 



/(u;, J, CO] ^ = /(a;, a;', CO - ^(a; - ^' - CO jd^fi^,^', f) 



■ (4.12) 



(4.13) 



Inserting 7j^^_|_i^_ into the definition of the three-particle distribution amplitudes one obtains the 
following expression 

03(0;', CO = 2(2 - D) {^a{oo',0 - *y(a;',CO) , (4-14) 
which is exactly the combination arising in the constraint derived from the equation of motion of the 
light quark in ref. [38]. At order ag, the other three-particle distribution amplitudes do not mix with 
</,_. 

One may use that to order the following relations hold [17]: 



duj 



dO-{uj; ji) 

dlogfi 
d(l)-{uj;n) ^ asin) 

d log /U Air 

and 



Slog fi 



0_(a;';/i)- J 



dlog n 



03(0;', e';/^), (4.15) 



J <koi}\LO,J-ix)4>.{J-ix) + j dw'd^'-i%{u,J,i'-ix)[^A-^v]{^\i'\lJ^) 



az_,3(.,.',e';M) _ .2e^,«(^,^',eo. 



d log /i ~ 47r 

taking into account that Z_^3 starts only at 0{as)- This leads to the anomalous dimension 7^1 3: 
.,(1) ,J d. .A - J^) 



(4.16) 

(1) . 



7i_^3(a;,a; ,C ;/x) = 7-;3('^>^^ 



CO' 

- Ca 



(Ca - 2Cf) 



-e(c -uj) + 



^'2 ^/ + _ ^ 



- m (®('^ - w ) - 6(0; - u; - O) 



(a;' + CO=^ i 

corresponding to the one-loop mixing between 0_ and ^a — *y. 



(4.17) 



J + 



5. Calculation in a general covariant gauge 



We have computed the mixing between gauge-invariant operators, and we could in principle have 
chosen any gauge to perform our determination of the renormalisation properties of the latter. We 
can check the validity of our previous computations by computing the nontrivial diagrams considered 
previously in a general covariant gauge, where we replace the Feynman-gauge gluon-propagator by 



ab 



g^iu - (1 - a) 



(5.1) 



A;2 

Our result should be gauge invariant, so that the parts proportional to (1 — a) should cancel. 
5.1 Two-parton external state and Z± 

First, we can repeat the computation of refs. [34,35], recalled in sec. |^, for the case of a two-parton 
external state. The gauge dependent part is: 

Mho± ^ -ig'sCpil -a) I ^v4^m^ {6{u; - k+ - /+) - 6iu; - k+)} , (5.2) 

j^v^^TT^u-^ {6{oj -k+ + /+) - 5{oj - k+)} , (5.3) 

-^^vi^TT''u-^5{oj-k+ + l+). (5.4) 

These integrals should be equipped with an infrared regulator (for instance, a gluon mass nig) in 
order to ensure that we keep only the ultraviolet divergences of interest here when we pick up the 
poles in e (otherwise, dimensional regularisation would treat both ultraviolet and infrared divergences 
of the integral as poles in e). 

In ref. [39], such integrals were considered with a particular focus on the integration over the 
different light-cone components d'^l —>■ ^dl+dl-(Pl±. Let us suppose that we want to integrate over 

There is a single pole, at I- = (m^ — ie + /^)//+, that we can always avoid by choosing the contour 
from above for /+ > and from below for /+ < 0. It seems to indicate that such integral should be 
0, which is incorrect. As proposed in ref. [39], a proper regularisation leads to the conclusion that 
an integration over the minus (plus) component results in a delta-distribution 5(/+) ((5(/_)), and one 
gets the equality: 

-Mio -k+± 1+) = 5{u; - k+) / (5.5) 



(2^)4/4-^^ ^ ^ J (27r)M4 

We see that Mho and Mlq vanish, whereas Mhl cancels the gauge-dependent part of the 
wave-function renormalisation for the heavy and the light quarks: 

ry 1 ry 1 

Z, = l + Ci.-i-[-l + (l-a)] Z;, = l + Ci.-^-[2 + (l-a)] (5.6) 
47r e 47r e 

Therefore, the gauge-dependent parts of the different contributions cancel and we have checked that 
the expression of Z± is indeed gauge independent. 

5.2 Three-parton external state and Z 3 

The issue becomes a little more involved if a three-particle state is considered. The complete formulae 
can be found in App. ^ but we can outline the pattern of cancellation for the gauge-dependent part 
among the various diagrams. 



We can identify the different gauge-dependent contributions in eq. ( 4.10| ). The diagrams (AAA), 
(i?44) and (C44), which vanish triviahy in the Feynman gauge, have to be taken into account, but 
their contributions can be shown to vanish through eq. (5.5). This is also the case for the diagrams 
(j414), (A24) and (A34). (-B13) and (-B14) remain finite as in the Feynman gauge. 

In analogy with the two-particle case, the diagrams (-B34), (i?35), (-B45) and (C34), (C35), 

1 /2 1 /2 

(C45) cancel with the gauge-dependent part of Zg' Z^' multiplied by B and C respectively. One 
has to pay attention to (C35) and (C45) that give additional contributions canceling each other. 
{B12), {B15), {B55) and (C12), (d5), (C55) cancel against the gauge-dependent part of Z^"^ Zg 
multiplied by B and C respectively. 

Finally, one is left with (A12), (A13), (A23), (523), (C13), (C14) and z]!'^ z]l'^ Z gZ^'^ x A. The 
sum of (C13) and (C14) is finite, and once eq. (|5.5|) is applied, only the following expressions remain: 





= -< — 


(^13) 


■ -JCa 

= [-J- 


(^23) 


■ Z^A 
= -^9s — 


(523) 


■ Z^A 

= -^9s — 



^^""^/ ■^Ji^{^i^-k+-^+)-^i^-''+)}^^±^T-ue+, (5.7) 

■^^^{S{u;-k+-q+)-6{u;-k+)}v7^^me+, (5.9) 

/d^l 1 1 
J^Ji— - ^+ - - - ^+)} ^4±TT''ue+. (5.10) 

Picking up the ultraviolet divergences from the integrals, we obtain finally for the sum: 

f% + Cf] (1 - a)-v4,TT^ue+— {6{uj - k+ - q+) - 6{uj - k+)} (5.11) 
which cancels exactly the (1 — a)-dependent part of the combination 

zl'^ZrZgZr X ^ = (l + £7 (^[1 + 2(1 - a)] - %[4 - (1 - a)])) x A (5.12) 

The calculation for is simpler, since the replacement of ^_ by r^j^ implies the absence of 

contributions proportional to v^r^j^r^_TT°'u. Therefore the additional terms of (C35) and (C45) 
vanish as well as the gauge-dependent contributions of the diagrams (C13) and (C14). Following the 
same lines in both cases, one can conclude that there is no gauge-dependence for the renormalisation 

of 6^ and (h_ . 



6. Conclusion 

In this paper, we have studied the mixing of both two-particle distribution amplitudes i?!>+ and 0_ 
with three-particle ones up to one-loop. Using the fact that RGE is a short-distance property of 
the operator, we used matrix elements of the operators with a quark-antiquark-gluon external state. 
Determining the ultraviolet divergences of the corresponding diagrams allowed us to recover the 
known one-loop self-mixing of (/>+ and , but also to determine the role of three-parton distribution 
amplitudes. We have established that (pj^ mixes only with itself, whereas (j)- does mix with [^^i — ^'y], 
and we have provided the corresponding anomalous dimension. Through the use of a general covariant 
gauge, we have checked that our results were indeed gauge invariant, providing further support to 
our expressions. 

We can relate our results to other comments on the 5-meson distribution amplitudes in the 
literature. For instance, the fact that <?!)_)_ does not mix with three-parton distribution amplitudes 



was already presented in ref. [40]. In this article, a computation similar to ours was sketched in 
the case of 4>^, with the conclusion (presented in eq. (2) of this reference) that the only ultraviolet 
one-loop divergence for 0+ is proportional to itself, whereas contributions proportional to higher- 
dimension operators have only infrared divergences. 

As mentioned in the introduction, the presence of 5{uj — uj')log{fi/uj) in the renormalization- 
matrices provides a radiative tail to (j)± falling off like (logw)/w for large oj. It requires one to 
consider either negative moments of the distribution amplitudes (/)_, (/>+, or positive moments with 
an ultraviolet cut-off [17,34,35,40,41]: 

(^^)±(^) = / do; a;^</>± (6.1) 



Jo 

On the contrary, it is interesting to notice the limit 

lim / dujuj^ z':^{{uj,uj',^') = N = 0, 1 (6.2) 

This is relevant for the calculation of the three-particle contributions to the moments: 

I duj uo^ = 1 + -r^ ( / duj' (j)-{uj') / duj uj^ z_ {uj,uj'; fi) (6.3) 

Jo 47r \J Jq 





Therefore, there is no contribution to the two lowest moments of 0_ from three-particle distribution 
amplitudes, which confirms the statement made after eq. (62) in ref. [35]. We have explicitly checked 
that this property of z^_^\ does not hold for higher positive moments (A^ > 2). 

In ref. [38] (see also refs. [42-44]) were derived two different relations between (/>_|_ and (/>_ on one 
hand and the four three parton-distribution amplitudes ^A, and on the other hand: 

uj(I)'_{lo) + (p+iuj) = I (to) (6.4) 
{uj-2A)(l)+{uj)+uj(l)^{uj) = J{uj) (6.5) 

/ (J) is an integro-differential expression involving — i^A + ^A and ^v)- The first relation 
comes from the equation of motion for the light quark and the latter one from the heavy quark (as 
suggested by the presence of the HQET parameter A = Mb — mi,). The use of the equation of 
motion of the heavy quark was criticised in ref. [33,35], because this is linked to the heavy-quark 
limit which does not commute with the light-cone limit. Moreover, this equation can be derived 
only if one leaves the light-cone limit, which is not needed for the first relation. One can notice that 
the shapes of the distribution amplitudes have been derived in the Wandura-Wilczek approximation, 
where three-parton distribution amplitudes are neglected (/ = J = 0), leading to rather unphysical 
shapes for the distribution amplitudes. If we assume that the three-particle distribution amplitudes 
mix separately into themselves [45,46], our work shows that the renormalisation-scale dependence 



of eq. (6^), derived from the light-quark equation of motion, is satisfying: by applying d/dlogfj, 
to the equation of motion, both (p^ and — yield a term proportional to — '^v- The 
contributions proportional to the two-particle distribution amplitudes were shown to cancel in the 
Wandzura-Wilczek approximation in App. D of ref. [35] (eq. ( |6.4D was also shown to hold in a specific 
non-relativistic model for (/)+,(/>_, "^y in App. C of the same reference). Eq. (6.5) does not seem to 



have such a satisfactory renormalisation-scale dependence, which would add to the various criticisms 
raised against this equation (see ref. [45] for further discussion of this issue). 



More generally, the influence of three-particle distribution-amplitudes on (j)^{uj;^) requires one 
to model them. However, the only available models [31] assume ^^(cj,^) = ^!v{'-^-X) ^-nd they yield 
no contribution to the evolution of ^/>-(a'; //). For practical calculations as well as for further model- 
building of distribution amplitudes beyond 0+, one needs the evolution kernel of the three-particle 
distribution amplitudes, which will be the subject of a future work [46]. 
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A. Extraction of poles related to UV divergences 

In this paper, we compute various integrals in dimensional regularisation to extract the e-poles 
related to UV divergences. The textbook procedure consists in a covariant anlaysis, where all space 
directions are treated on the same footing. Since we use light-cone coordinates with a privileged 
direction for the definition and discussion of the distribution amplitudes, we use a slightly less usual 
method which we apply on the illustrative integral: 



(27r)4'' ^^^^2 {i-kf 

where / is an arbitrary function of alone, corresponding to a gluon line (of momentum /) attached 
to a light-quark line (of incoming momentum k). Such an integral is needed already to compute the 
mixing of (})± into themselves {Zio)- We want to perform the integrals over and li_ in 4 — 2e 
dimensions and isolate the poles in e related to UV divergences. We therefore introduce a small mass 
m for the light quark to regularise (soft) IR divergences that are of not interest for the determination 
of the RG properties of the distribution amplitudes. 

We perform first the integral over Z_ by identifying the poles in the complex /_ plane : 

_ r di^di.dPi^ 1 1 

j 2(27r)4 •^^'+^Z2 + iO+Z2-2A;-/ + zO+ ^ ' 

which are = — ■iO"'")/l_|_ and I- = k- + {{l±_ — k±))'^ — zO"'")/(Z+ — whose positions with 
respect to the real axis depend on the value of l+. If 1+ is negative or larger than k+, the two poles 
sit on the same side, and thus the contour integral yields 0. If < Z+ < the two poles are on 
different sides and one gets a non-vanishing contribution, for instance by closing the contour in the 
lower half-plane and thus picking up the first pole (associated with the /^-denominator): 

Then one can perform the integral over the 2 — 2e transverse dimensions, which yields the result: 



2 ;2 



m I 



(A.4) 



This expression yields a single pole in e corresponding to the UV-divergent part of the integral, 
which enters Zlq- The same procedure is applied to all the diagrams, with sometimes more involved 
integrals (up to four propagators), leading to the results quoted in the present article. 



B. Gauge dependence of the diagrams in a general covairiant gauge 



In this appendix, we collect the integrals from the different diagrams that are proportional to the 
gauge parameter (1 — a). For the well known cases (B12), (B15), (B55), (C12), (C15) and (C55), 
corresponding to the vertex renormalisation, the integrals have already been carried out. 
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a 1 1 
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As mentioned at the end of Sec. 6, the terms proportional to ^i^^i^j. vanish for (f)+. In addition, all 
the integrals must be understood with an infrared regulator, since we are only interested in their 
ultraviolet behaviour. 
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